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^ . Abstract 

a^ 

We calculate the masses of the lowest lying eigenstates of improved SU(2) and 
, SU(3) lattice gauge theory in 2+1 dimensions using an analytic variational approach. 

The ground state is approximated by a one plaquette trial state and mass gaps are 
, calculated in the symmetric and antisymmetric sectors by minimising over a suitable 

' basis of rectangular states. 



1 Introduction 

^ I Lattice gauge theory (LGT) [l| is a non-perturbative technique used to regulate the ul- 

traviolet divergences of non-abelian gauge theories. There are two standard formulations 
of LGT: the Lagrangian formulation and the Hamiltonian formulation. 

X 

• To date most work on LGT has been carried out in the Lagrangian formulation for which 

Monte Carlo techniques are readily applicable. This approach involves replacing infinite 
and continuous space and time by a finite and discrete lattice of points. The Hamiltonian 
formulation developed by Kogut and Susskind [^] differs by leaving time continuous. In 
this framework the problem of gauge theory reduces to a many body problem and one is 
interested in calculating ground state and excited state configurations of the field defined 
on the lattice. The relation to continuum physics is realised in the limit in which the 
lattice spacing vanishes and the extent of the lattice becomes infinite. 

Many body techniques, in particular plaquette expansion ^, [s], ^] and the coupled clus- 
ter method 0, |8|, ^ |lO|, have been applied to Hamiltonian LGT for a number of years. 
Not only do such calculations serve as a check on the more common Monte Carlo tech- 
niques of the Lagrangian approach, they also provide more immediate access to excited 
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states |jTT|. In addition, the problem of finite density QCD may be more easily handled in 
the Hamiltonian approach |12|. 



In recent years an improvement program has been undertaken in the Lagrangian approach. 
The aim has been to systematically correct the discretisation errors arising in the lattice 
Lagrangian and hence decrease the computational cost of lattice calculations [13|. More 



recently this improvement program has moved to the Hamiltonian formulation |14]. In 
a recent paper we demonstrated how improved Hamiltonians can be constructed [^]. 
In this paper we explore the use of improved Hamiltonians in variational calculations of 
SU(2) and SU(3) vacuum energy densities, specific heats and exited state energies on a 
2+1 dimensional lattice. A key development in this paper is the application of analytic 
techniques to SU(3) LGT. 

The outline of this paper is as follows. In section 2 we introduce a generating function 
for SU(3) which allows us perform analytic calculations of energy densities and specific 
heats. In section 3 we calculate the energies of the lowest lying eigenstates of SU(2) and 
SU(3) pure gauge theory in 2+1 dimensions and compare the results of unimproved and 
improved calculations. Section 4 contains our conclusions. 

2 Vacuum Energy Density and Specific Heat 
2.1 Introduction 

In this section we calculate vacuum energy densities and specific heats for the Kogut 
Susskind, improved and tadpole improved Hamiltonians. Improved Hamiltonians are eas- 
ily constructed by adding appropriately weighted gauge invariant terms to the Kogut 



Susskind Hamiltonian [15|. For SU(A^) pure gauge theory with coupling on a lattice 
with spacing a, we define a class of Hamiltonians in terms of the link variables Uj{x) and 
electric fields Ei{x) (which are themselves defined in the appendix), and two parameters 
K and no which define the particular Hamiltonian in question: 



H{k,uo) 



a ^ 



(1 - K)Ei{xf + -^Ei{x)Ui{x)E^{x + ai)Ul{x) 



+^ E {il+^'^)nj{x)-'^[Ri,{x) + R,i{x)] 

x,i<j 



where 



P^jix) 



^ReTi \Ui{x)Ujix + ia)U^^{x+ja)uUx] 
Uq I- ■' 

(D) 



Nut 



and 



Rij{x) 



Re 



(3) 
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are the plaquette and rectangle operators respectively. The Kogut-Susskind and 0(0^) 
classically improved Hamiltonians are given by Hks = -^(0,1) and Hj = H (1/6,1) re- 
spectively. The tadpole improved Hamiltonian is given by Hxi = H(1/6,uq), where uq is 
now the mean link, defined by 

and evaluated self-consistently, as a function of g'^. Details of the mean link calculation 
are given in section |2.3| . 



We approximate the ground, or (perturbed) vacuum, state, with energy Eq, with the one 
plaquette trial state, first suggested by Greensite pi]. 



exp 



Zp + Zp) 



[0). 



(5) 



Here, |0) is the strong coupling vacuum defined by -Ef (a;)|0) = for all i, x and a 
1,2,... ,N'^ -1. The p-sum extends over all plaquettes on the lattice and 



Zn 



Tr 



Ui{x)Uj{x + ia)Uj{x + ja)uUx 



(6) 



denotes the plaquette labelled by p which joins the lattice sites x, x + ia, x + {i + j)a 
and X + ja. Zp is the trace of the link operators in the opposite direction. For a given 
coupling the variational parameter c is fixed by minimising the vacuum energy density, 
which, after some algebra making use of eqn(|47|) and eqns(55 — 59) from the appendix, can 
be expressed in terms of the expectation values of plaquettes and rectangles as follows: 



eo 



iV2 



2P 



l\ 2/3(1 + 4k) 
c - 




Nut 



(7) 



Here Np is the number of plaquettes on the lattice and f3 = N/g is the inverse coupling. 



2.2 Analytic Techniques for SU(2) and SU(3) 

The expression for the SU(A^) vaccuum energy density given by eqn(0) allows us, in prin- 
ciple, to fix the variational parameter c for a given f3 in any number of dimensions. The 
difficulty lies in the calculation of the plaquette and rectangle expectation values. Monte 
Carlo simulations can be used to calculate these expectation values in any number of 
dimensions. However, for the special case of 2+1 dimensions, the calculation can be car- 
ried out analytically. This is because the change of variables from links to plaquettes 
has unit Jacobian |17]. Consequently, the plaquettes on the two dimensional lattice are 
independent variables, which leads to: 




lUidUiUpe' 



c{Zp+Zp) 



(8) 
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and 



(9) 



Eqn(^) makes use of the independence of plaquettes and eqn(p^. For the case of SU(2), 
analytic expressions for the plaquette expectation value in terms of modified Bessel func- 
tions (/„) have been used in variational calculations for the last 20 years: 



a 



2/2 (4c) 
/i(4c) 



(10) 



We derived the corresponding SU(3) result, which has not been used, to our knowledge, 
in the Hamiltonian approach, and then discovered that it follows simply from a paper of 
Eriksson et al [^], in which the SU(3) integral, 



dUe 



Tr((7Mt+;7tM) 



(11) 



is calculated for arbitrary 3x3 matrices M. Following their analysis, but treating M and 
as independent variables, with M = cl and = dl, leads to an expression for the 
SU(3) generating function: 



yic,d) 



— f dz-^— 
n Jr ziz 



cd) 



3/2 



(12) 



Here Ji is the first Bessel function and F is a closed contour in the complex plane including 
the pole at z = but excluding the pole aX z = cd. To evaluate this contour integral we 
expand the integrand in power series about the pole at z = and use Cauchy's integral 
theorem to eventually obtain. 



LXJ -| ft, 



3fc + 3 
k-l 



■cd) 



k-L 



c' + d 



(13) 



This generating functional is extremely useful for Hamiltonian lattice calculations^ In 
principle, it allows us to investigate 2+1 dimensional pure SU(3) gauge theory analytically. 
The calculation of various matrix elements for all couplings, reduces to a exercise in 
differentiation. For example: 



(Zp) 
(ZpZp) 



a)- 

1 a™+"3; 



'\_dy_ 
y dc 

"i d'^y 

ydcdd 



d=c 



, and 



- d=c 



(14) 



y dc'^dd" 



^The SU(A'^) generalisation of this result, which we have recently obtained, will be reported in a future 
publication. 
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Figure 1: Analytic and Monte Carlo calculations of the SU(3) vacuum energy density 
in 2+1 dimensions for /3 = 4.0. It should be noted that the errors in the Monte Carlo 
calculation are smaller than the plot points. 



As a check on this series we calculate the strong coupling limit (c = 0) of {Z^Zp). By 
differentiating eqn(13), we observe that the only non-zero strong coupling matrix elements 



occur when n + 2m = mod 3. For this case we have, 

IZmzn) _ = V f ^^^^ ^ (l^) 

Here the sum runs over all integers B±!li < k < min "+^^'" ^ m+2n ^ ^ Thig strong coupling 
result has an equivalent combinatoric interpretation as the number of times the singlet 
representation appears in the direct product of m 3 and n 3 representations A simple 
Mathematica code verifies the agreement of the two approaches. 



A demonstration of the accuracy of the series away from the strong coupling limit is 
provided by a Monte Carlo calculation of the vacuum energy density at various couplings. 
The example of /? = 4.0 is shown in fig |l] and indicates close agreement between the Monte 
Carlo and analytic calculations. 



2.3 Results 

To calculate the energy density and variational parameter as functions of f3 we proceed as 
follows. For the Kogut-Susskind and classically improved cases, we simply minimise eo for 
a given value of (3. The value of c at which eo takes its minimum defines c as a function of 



The tadpole improved calculation is more complicated. This is because the mean pla- 
quette, from which uq is calculated, depends on the variational state {(po). To complicate 
matters, the variational state in turn depends on the energy density and hence uq. Such 
interdependence suggests the use of the following iterative approach. For a given /? and 
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Figure 2: Analytic calculation of the unimproved, improved and tadpole improved vacuum 
energy density in 2+1 dimensions for SU(2) and SU(3). 



initial approximation of uq we minimise the energy density to fix the variational state |i;^o)- 
Using this trial state, we update uq for use in the next iteration. This process is iterated 
until convergence is achieved, typically after only a few iterations. 

The results of the Kogut-Susskind, classically improved and tadpole improved SU(2) and 
SU(3) vacuum energy density calculations are shown in fig |^. The corresponding varia- 
tional parameters c(/3) are shown in fig ^. The correct strong and weak coupling behavior 
is observed in each case. The differing gradients in the weak coupling limit highlight the 
fact that when using an improved Hamiltonian one is using a different renormalisation 
scheme to the unimproved case. 

In practice the SU(3) generating function is truncated. The dependence of the variational 
parameter on various truncations of the fe-sum in eqn(13) is shown in fig ^. We see that 
convergence is achieved up to /? ~ 13 when keeping 20 terms. Further calculations show 
that when keeping 50 terms convergence up to /3 « 30 is achieved. 



In addition to the vacuum energy density we can also calculate the specific heat: 



Cv 



dp' 



•2 ■ 



(16) 



The results are shown in fig ^. We note that the location of the peak indicates the region 
of transition from strong to weak coupling. For an improved calculation one would expect 
the peak to be located at a smaller /5 (corresponding to a larger coupling) than for the 
equivalent unimproved calculation. We see that this is indeed the case for both SU(2) 
and SU(3), with the tadpole improved Hamiltonian demonstrating the largest degree of 
improvement. 
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(a) SU(2) (b) SU(3) 



Figure 3: Analytic calculation of the unimproved, improved and tadpole improved varia- 
tional parameter in 2+1 dimensions for SU(2) and SU(3). 




5 10 15 20 

Figure 4: Analytic calculation of the unimproved SU(3) variational parameter in 2+1 
dimensions, truncating the /c-sum of y{c, d) at k = kmax- 
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Figure 5: The unimproved, improved and tadpole improved specific heat in 2+1 dimensions 
for SU(2) and SU(3). 



3 Mass Gaps 

3.1 Introduction 

Having fixed the chosen trial vacuum \4>q), in this section we turn to investigating excited 
states. Our aim is to calculate the lowest lying energy eigenstates of the Hamiltonians 
described by eqn(|) for both SU(2) and SU(3). 

We follow Arisue [^] and expand the excited state in the basis consisting of all 
rectangular Wilson loops {10}^'"' that fit in a given square whose side length L^ax 
defines the order of the calculation. In order to ensure the orthogonality of |(/>o) and 
we parameterise the excited state as follows 

r r2 

J-'max max 

= J2 ■■^i\^^^)= E (17) 

n,m=l 1=1 

with, 

iO = J2iWi{^)-iWi{x))]\^o). (18) 

Here {Wi{x)) is the expectation value of Wi{x) with respect to the ground state |0o) and 
the convenient label I = {n — l)Lmax + has been defined to label the nx m rectangular 
state. We define the particular form of Wi{x) to reflect the symmetry of the sector we 
wish to consider. For SU(3) we take Wi{x) = Tv[wi{x) it wj{x)] for the symmetric sector 
and antisymmetric sector. To avoid over-decorated equations, the particular Wi{x) in use 
is to be deduced from the context. Here wi{x) is the rectangular Wilson loop joining the 
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lattice sites x, x + nai, x + nai + maj and x + naj, with 



n 



l-l 



+ 1 and m = I — Lr, 



l-l 



(19) 



Here [k] denotes the integer part oi k. In order to calculate excited state energies we 
minimise the mass gap (the difference between the excited state and ground state energy) 
over the basis defined by a particular order Lmax- To do this we again follow Arisue [^ ] 
and define the matrices 



Nu 



1 



{l\H-Eo\l') 



(20) 



and 



Dw = ^m = j2[{wi{x)Wi,m-{Wi{x)r{Wi,m]. (21) 

P X 

Extending the calculation to the general improved Hamiltonian H{k, uq) and making use 
of eqns(|56| — 59) from the appendix gives 



Nn 



'2a 



EE (1 - ^){[En^),wi{x')] [Er{x),w,m 

i,x x' K 



(22) 



+4 ( iKix), W^ix')] \e^{x + ai), Wi>{0) 

^0 

To minimise the mass gap over a basis of a given order we make use of following diago- 

(23) 



nalisation technique |21]. We first diagonalise D, with 

S^DS = V^, 

where V is diagonal. The n-th lowest eigenvalue of the modified Hamiltonian 

H' = VS^NSV, (24) 

then gives an upper bound for the mass gap corresponding to the n-th lowest eigenvalue 
of the Hamiltonian, Am.„. 



3.2 Calculating Matrix Elements 

Having described the minimisation process we now detail the calculation of the matrix 
elements Nu' and Du' . Our aim is to reduce Nu' and Dn' to polynomials of one plaquette 
matrix elements. This, again has been done for the case of SU(2) by Arisue [20|. Here we 



retrace his calculations for the case of SU(3) and extend them to incorporate improved 
Hamiltonians. We start with Du'. 

Taking elementary plaquettes as our independent variables, it is easy to show that the only 
non-zero contributions to Du/ occur when the rectangles I and /' overlap. As an example 
of a contribution to Dni , consider ADu' ; the case where Ni^i' plaquettes are contained by 
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both rectangles (these are the overlap plaquettes) and Ni plaquettes are contained by the 
rectangle /. In order to calculate such matrix elements we rely heavily on the orthogonality 
properties of SU(A^) characters. We are interested in calculating SU(A^) integrals of the 
form 

J dUpe^'-^^hriUpV), (25) 

where Up is a SU(A^) plaquette variable and F is a product of any number of plaquettes 
not including Up. Here Xr{U) denotes the character corresponding to the representation 
r. For SU(2), r = 0, 1/2, 3/2, . . . and for SU(3), r = (A, ^) where A denotes the number of 
boxes in the first row of the Young's Tableaux describing the representation and /i is the 
number of boxes in the second row. Similarly, for SU(A^), r = (ri,r2, . . . ,r7v-i)- 



Performing a character expansion of the exponent in eq. ^ gives: 

dUpe^^'^^'hriUpV) = T.J dUpCr^Xr'{Up)Xr{UpV). 



(26) 



This is simply a generalisation of a Fourier expansion. Here, the coefficient c^/ is given by: 



dUpXr'{Up)e^^^^\ 



Now using the orthogonality relation, 

dUpXr'{UpV)Xr(Up) 



d, 



■5r'rXr(y), 



where dr is the dimension of the representation r, we obtain: 

dUpe^^^^^XriUpV) = ^Xr{V) J dUpXr{Up)e'^^^-\ 



(27) 



(28) 



(29) 



This result allows us to integrate out a single plaquette from an extended Wilson loop in a 
given representation r. To complete the calculation we need to relate characters to traces 
of group elements. This can be done using Weyl's character formula |22|. For SU(3) the 
dimensions and characters corresponding to the first few representations are given by: 



XioiU) 
Xn{U) 

X2o{U) 

X2l{U) 



TrC/ 
1 



{TiUy - Tr(C/' 



TrU^ 



{TiUf + Tr(C/2)l = {TiUf - TiU^ 



TtUTtU^ - 1 



dio 
dii 

d20 

d2i 



(30) 



{TrUf -Tt{U^) 

Here we have made use of the SU(3) identity 

Tr(C/2) = {TrUf - 2TiUl (31) 

Applying eqns(|2^ — |2l|) allows us to obtain analytic expressions for each contribution to 
Dill. For the case of AD/// described earlier, we have 



AA/' = ^Fz{Ni + Ni, - 2Nini') [Fz2{Nini') + FzziNmi')] - 4Fz{Ni)Fz{Ni,), (32) 
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where the character integrals are given by: 



Fzin) 
Fz2{n) 



n-l 



n-l 



(Zp — Zp)"' 



n-l 



{Zp}"^ , and 



1+^ 



n—l 



{{ZpZp)-lf^ 



(33) 
(34) 
(35) 



We now move on to the calculation of Nui. It is easy to show that the only non-zero 
contributions occur when there is at least one common link and an overlap between the 
rectangles. The improvement term (the second term in eqn(p2|)) only contributes when 
the two rectangles share at least two links in a given direction. Consider the contribution 
ANiii to Nil' in which there are Li common links and L2 common strings of two links in 
a given direction. Again we suppose Ni plaquettes are enclosed by rectangle / and that 
there are Ni^i' common plaquettes. Making use of eqn (p9|) and eqns(56 — 59) from the 
appendix we obtain 



AiV„, = ^^FziNi+Ni,-2Nini' 



;Fz2{Nini')-2Fz{Nini')- 



-^+lFzz{Nini' 



with 



L 



k)Li H 2-^2- 



(36) 



(37) 



Having calculated individual contributions to Du/ and Nw, the completion of the calcu- 
lation requires counting the possible overlaps of a given type. The approach described in 
this section is easily extended to SU(A^). However, analytic expressions for the character 
integrals are not available for this case. 

3.3 Results 



For each SU(3) calculation we have kept 80 terms in the /c-sum of eqn(13) giving conver- 
gence up to /5 = 50. The generation of A'^;;' and Dni and implementation of the minimisa- 
tion process was accomplished with a Mathematica code. 

For the case of 2-1-1 dimensions we expect Am/g'^ to become constant in the scaling region. 
The convergence of the massgaps with Lmax for the cases of SU(2) and SU(3) is illustrated 
in figs 5(a) and 6(b). We notice that for the case of SU(3) only small improvements to 



scaling are gained by extending the calculation beyond order 8. This suggests that a more 
complicated basis (including, for example, loops covered more than once) is required to 
simulate SU(3) excited states than for the case of SU(2). 

From fig |^, we see that Amf / is approximated well by a constant in very large scaling 
regions for the lowest lying eigenstates for both SU(2) and SU(3). The scaling behavior 
becomes significantly worse for the antisymmetric sector and for higher energy eigenvalues. 
This is because the simplistic form of our excited state wave function is not sufficient to 
reproduce the plaquette correlations required to simulate these higher order states. One 
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Figure 6: The unimproved 2+1 dimensional symmetric massgaps for SU(2) and SU(3). 

would expect the simulation of higher order eigenstates to improve by including more 
complicated loops in our expansion basis or by using a more complicated ground state. 
The continuum limit results of excited states for SU(2) and SU(3) results are given in 
tables |-^. 

For the unimproved SU(2) case, the lowest two eigenstates agree precisely with the cal- 
culation of Arisue . This serves as a check on our counting in calculating the possible 
overlaps of excited states. Our calculation is in disagreement with that of Arisue at the 
3rd eigenstate. It appears that our 4th eigenstate corresponds to Arisue's 3rd and that 
our third eigenstate does not appear in his results. The reasons for this are not clear. 



The results for the SU(3) symmetric massgap (in units of g^) are to be compared to calcu- 



lations by Luo and Chen 2.15±0.06 [||, Samuel 1.84±0.46 |3l and Teper 2.40±0.02 [||]. 
It should be stressed that although our result of 3.26520 it 0.00009 is considerably higher, 
it provides a precise upper upper-bound on the exact result. By including more compli- 
cated loops in the expansion basis one would expect to reduce this upper bound. This is 
emphasised by the fact that when using only square basis states the result is considerably 
higher. 



When compared to equivalent unimproved calculations, the improved and tadpole im- 
proved massgaps approach scaling faster as (5 is increased. This is evident in figs ^ and |^ 
and is expected since, for an improved calculation one is closer to the continuum limit 
when working at a given coupling. However, for each improved calculation the scaling 
behaviour is not as precise as the equivalent unimproved calculation. This is clear from 
tables H — 13[ A possible reason for this is that the one plaquette trial state used here does 
not allow for direct contributions from the improvement term in the kinetic Hamiltonian. 
For this term to contribute directly one would need a trial state which includes Wilson 
loops extending at least two links in at least one direction. 
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Figure 7: 2+1 dimensional massgaps for SU(2) and SU(3) with L^nax = 25. 
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Unimproved 


Improved 


Tadpole Improved 


Amf 


2.05691 ± 0.00002 


2.0897 ± 0.0003 


2.0965 ± 0.0006 


Amf 


3.645 ± 0.001 


3.685 ± 0.001 


3.6953 ± 0.0009 


Amf 


4.5202 ± 0.0004 


4.574 ± 0.004 


4.583 ± 0.004 


Amf 


5.133 ±0.003 


5.177 ±0.004 


5.189 ± 0.004 


Amf 


5.867 ±0.006 


5.932 ± 0.008 


5.943 ± 0.008 



Table 1: Upper bounds on the lowest lying SU(2) massgaps (in units of g^) computed 
with various Hamiltonians in 2±1 dimensions. The unimproved, improved and tadpole 
results are calculated in the respective scaling regions 13.5 < /? < 30.0, 9.9 < /? < 30.0 
and 9.25 < /? < 30.0. 





Unimproved 


Improved 


Tadpole Improved 


Amf 


3.26506 ± 0.00007 


3.3201 ± 0.0003 


3.3350 ± 0.0008 


Amf 


6.194 ±0.003 


6.273 ± 0.003 


6.303 ± 0.002 


Amf 


7.396 ± 0.010 


7.536 ± 0.008 


7.575 ± 0.006 


Amf 


8.844 ± 0.007 


8.936 ± 0.008 


8.980 ± 0.006 


Amf 


9.87 ±0.02 


10.04 ±0.01 


10.10 ±0.01 



Table 2: Upper bounds on the lowest lying symmetric SU(3) massgaps (in units of g"^) 
computed with various Hamiltonians in 2±1 dimensions. The results are calculated in the 
scaling region 18.5 < (3 < 50.0. 



The improved SU(2) massgap can be compared to the coupled cluster calculation of Li et 
al |2^. Their result (in units of g"^), Amf = 1.59, is again significantly lower than our 
result 2.0897 ± 0.0003. The difference is attributable to the different choices of ground 
state used. While our calculation makes use of the simple one plaquette ground state, the 
coupled cluster calculation of Li et al uses a more accurate ground state wave function 
consisting of an exponential of a sum of extended loops. 

We can also calculate massgap ratios. For example the antisymmetric-symmetric massgap 
ratio for SU(3) is given by 



Amf 
Amf 



1.6173 ± 0.0008 Unimproved 

1.6177 ± 0.0006 Improved (38) 
1.6182 ± 0.0005 Tadpole improved. 



These results lie between the unimproved results of Chen et al, 1.6989 [27] and Teper, 
1.50 ±0.03 pH]. 
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Unimproved 


Improved 


Tadpole Improved 


Amf 


5.281 ±0.003 


5.371 ± 0.002 


5.396 ± 0.001 


Ami 


7.831 ± 0.007 


7.947 ± 0.006 


7.986 ± 0.004 


Ami 


7.40 ± 0.01 


8.87 ±0.01 


8.922 ± 0.009 


Ami 


10.32 ±0.01 


10.46 ±0.01 


10.514 ±0.009 


Ami 


11.06 ±0.02 


11.27 ±0.02 


11.33 ±0.01 



Table 3: Upper bounds on the lowest lying antisymmetric SU(3) massgaps (in units of g'^) 
computed with various Hamiltonians in 2±1 dimensions. The results are calculated in the 
scaling region 18.5</3<50.0. 



4 Conclusion 

In this paper we have extended the analytic techniques of 2±1 dimensional Hamiltonian 
LGT, traditionally used for SU(2), to the case of SU(3). Impressive scaling is achieved over 
an extremely wide range of couplings for the lowest energy eigenstates in the symmetric 
and antisymmetric sectors. Our results use a one plaquette trial state and a basis of 
rectangular states over which excited state energies are minimised. Such choices allow us 
to use analytic techniques for both SU(2) and SU(3). Our results provide upper bounds on 
SU(2) and SU(3) unimproved, improved and tadpole improved glueball masses all of which 
are above current estimates. We expect that these upper bounds will lower with the use of 
more complicated states in the simulation of both the ground state and the excited state. 
For such choices of states the analytic techniques used here are not as readily applicable. 
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A Commutation Relations 

In improved Hamiltonian LGT calculations, one encounters matrix elements of the form: 

(0olETr[^»(a=)-E,(x)]|(/>o) and {(|)o\J2^^^[Ei{x)U^ix)E^ix + ai)U^^ix)]\cPo). (39) 

The first of these is easily handled. One simply writes the electric field operators Ei{x) in 
terms of their components with respect to the Gell-Mann matrices A", a = 1, . . . , N'^ — 1, 

Ei{x) = E^{x)X''. (40) 

Making use of the SU(Ar) relation 

TV(A"A'^) = ^6a/3, (41) 



we have, 

{ct>o\j2'^[Eii^)Eii^)]\^o) = ^j2i^o\Enx)Er{x)\4>o). m 



x.t 
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Let the trial state |i?!>o) have the form |0o) = C^lO), where 5 is a function of hnk variables 
and 5"^ = S. Then 



(43) 



Making use of the Baker-Hausdorff formula we can derive the following result, 



e^E^{x)e- 



Ef{x) - [Ef{x),S] + -[[Enx)^SlS] + ... 
E'-{x)-[E2{x),S]. 



(44) 



The last line follows from the fact that S" is a function of link variables and the commutation 
relation: 



[Et{x),Uj{y)] = 5,,6^yX''Ui{x). 
Some rearrangement of eqn(|4^) leads to the useful result, 

'Er{x),e'] = [Enx),S] = [Enx),S] , 
which we apply in eqn(|43|) to obtain, 

Y,{<Po\Enx)Et{xMo) = -j2m[Enx),s][Er{x),s]\^o) 

= lY.(M[Er{x\[Erix),s]]\c^o). 



(45) 



(46) 



(47) 



In the last line we have used integration by parts [28|. 

The second of the matrix elements in eqn(^) is more difficult to handle. For this case it 
is convenient to expand Ei{x + ai) in the basis Xf{x + ai): 



Ei{x + ai) = ~\1{x + ai)Ef{x + ai), 



where. 



\f{x + ai) = ul{x)X^Ui{x) 



(48) 



(49) 



We now derive the commutation relations between Ef[x) and Uj{y). We start with the 
analogous relation to eqn(^T|), 



Tr 



A?(^)Af(rr) 



(50) 



which follows trivially from eqn(49). This relation allows us to invert eqn(pq) and write 
E" in terms of E"' 



ETix) 



2Tr 
2Tr 



Ei{x)Xf{x) 
X?{x)XP] E^{x). 



(51) 



16 



We can form a similar relation between A" and A" as follows. Since Af (a;) € SU(A^), we 
can expand in the Gell-Mann basis as follows, 



(52) 



Here c^^ {x) are constants which we now determine. Multiplying eqn(52) throughout by 
A^, tracing and making use of eqn(pT|) gives 



cfix) = 2Tr 



(53) 



From eqn(p]) we see that Ef(x) = c^^ {x)E^ (x). The commutation relations follow 
immediately, 

'Er{x),U,{y)] = cf{x)[E^{x),U,{y)]=cf{x)X('5,,6^yU,{x) 

= Sijd^yXUxmx). (54) 

These commutation relations enable the simplification the second of the matrix elements 



in eqn(|39|) using the basis defined in eqn([48D. Making use of eqns ( |5l| ) and (|54D leads to 
the following simplification: 

(</>o|$^Tr [£;i(a;)C/i(a;)i5;i(a;+ai)?7;(a=)] |(/.o) = \ Y.{M [^^"(a;), l</'o).(55) 
In order to calculate such matrix elements the following results are useful 



E^{x), Ef{x + ai),Ui{x)U^{x + ai) 



2N 



'-Ui{x)Ui{x + ai), (56) 



[E^ix), {U,ix)}AB] Kix + ai), {U,{x + ai)}cD 
I 
2 

E^{x), {u}{x)}ab\ + ai), {U^{x + ai)}cD 



^6Bc{Ui{x)}AB'{Ui{x + ai)}B'D - j^{Ui{x)}AB{Ui{x + ai)}cD, (57) 



-l{uhx)}cBmx + ai)}AD + ^{uhx)}ABmx + ai)}cD, (58) 
[Efix), {Ui{x)}AB] [Efix + ai), {Ulix + ai)}cD 

-lmx)}cB{Uj{x + ai)}AD + ^mx)}AB{ul{x + ai)}cD. (59) 

Here {X}^^ denote the colour indices of the matrix X and an implicit sum over repeated 
colour indices is understood. These results follow simply from the commutation relations 
of eqns (^) and (^^, and the following SU(A^) formula: 



Q a 1 1 
^AB^CD = 2^AdSbC - l^^AB^CD- 
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